10.2 KB · ⇊ 3 parts · ⤓ Micron 10.2 KB · ⤓ HTML 73.4 KB

top

Mit Hilfe der gaußschen Trapezformel (nach Carl Friedrich Gauß) ist es möglich, die Fläche eines einfachen Polygons zu berechnen.cite-ref-1[1] Dabei wird jeder Polygonkante ein Trapez (siehe Bild) zugeordnet, dessen Flächeninhalt sowohl positiv als auch negativ sein kann. Negative Flächenteile kompensieren außerhalb des Polygons liegende Teile positiver Trapeze.

Eine Variation der Trapezformel ist die Dreiecksform, deren Analogon für stückweise glatte Kurven die Sektorformel von Leibniz ist.

• Das Prinzip und die Formel

• Beispiel und Schnürsenkel-Schema

• Einzelnachweise

• Weblink

• Literatur

Es wird vorausgesetzt, dass die Punkte P i = ( x i , y i ) , i = 1 , . . . n {\displaystyle P_{i}=(x_{i},y_{i}),i=1,...n} des Polygons im ersten Quadranten eines kartesischen Koordinatensystems liegen und die Punkte mit wachsendem i {\displaystyle i} im mathematisch positiven Sinn (Gegenuhrzeigersinn) durchlaufen werden. Aus praktischen Gründen wird angenommen, dass P n + 1 = P 1 {\displaystyle P_{n+1}=P_{1}} ist. Der Kante P i P i + 1 {\displaystyle P_{i}P_{i+1}} wird dann der Flächeninhalt

A i = 1 2 ( y i + y i + 1 ) ( x i − − x i + 1 ) {\displaystyle A_{i}={\tfrac {1}{2}}(y_{i}+y_{i+1})(x_{i}-x_{i+1})}

des Trapezes ( x i , y i ) , ( x i + 1 , y i + 1 ) , ( x i , 0 ) , ( x i + 1 , 0 ) {\displaystyle (x_{i},y_{i}),(x_{i+1},y_{i+1}),(x_{i},0),(x_{i+1},0)} zugeordnet. Ist x i < x i + 1 {\displaystyle x_{i}<x_{i+1}} so ist A i {\displaystyle A_{i}} negativ, im anderen Fall positiv oder A i = 0 {\displaystyle A_{i}=0} falls x i = x i + 1 {\displaystyle x_{i}=x_{i+1}} ist. In der Zeichnung ist die Orientierung der Kanten durch Pfeile gekennzeichnet. An der Farbe der Pfeile ist das Vorzeichen der jeweiligen Trapezfläche zu erkennen: rot steht für A i < 0 {\displaystyle A_{i}<0} , grün für A i > 0 {\displaystyle A_{i}>0} . Im ersten Fall heißt das Trapez negatives Trapez, im zweiten Fall positives Trapez. Die negativen Trapeze löschen die außerhalb des Polygons liegenden Flächenteile positiver Trapeze. Am einfachsten ist dies an dem Beispiel eines konvexen Polygons (im Bild oben) zu erkennen: Der Flächeninhalt des Polygons ist gleich der Summe der Flächeninhalte aller positiven Trapeze (mit grünen Kanten) minus den Flächeninhalten aller negativen Trapeze (mit roten Kanten).

Für den Flächeninhalt, des von dem Polygon P 1 , . . . , P n {\displaystyle P_{1},...,P_{n}} eingeschlossenen Gebiets ergibt sich also

A = 1 2 ∑ ∑ i = 1 n ( y i + y i + 1 ) ( x i − − x i + 1 ) {\displaystyle A={\frac {1}{2}}\sum _{i=1}^{n}(y_{i}+y_{i+1})(x_{i}-x_{i+1})}

= 1 2 ( ( y 1 + y 2 ) ( x 1 − − x 2 ) + ⋯ ⋯ + ( y n + y 1 ) ( x n − − x 1 ) ) {\displaystyle \quad ={\frac {1}{2}}{\Big (}(y_{1}+y_{2})(x_{1}-x_{2})+\cdots +(y_{n}+y_{\color {red}1})(x_{n}-x_{\color {red}1}){\Big )}}

Multipliziert man die Klammern aus und beachtet ∑ ∑ i = 1 n x i y i = ∑ ∑ i = 1 n x i + 1 y i + 1 {\displaystyle \sum _{i=1}^{n}x_{i}y_{i}=\sum _{i=1}^{n}x_{i+1}y_{i+1}} , erhält man die Determinantenform der Flächenformel:

A = 1 2 ∑ ∑ i = 1 n ( x i y i + 1 − − x i + 1 y i ) = 1 2 ∑ ∑ i = 1 n | x i x i + 1 y i y i + 1 | {\displaystyle A={\frac {1}{2}}\sum _{i=1}^{n}(x_{i}y_{i+1}-x_{i+1}y_{i})={\frac {1}{2}}\sum _{i=1}^{n}{\begin{vmatrix}x_{i}&x_{i+1}\\y_{i}&y_{i+1}\end{vmatrix}}}

= 1 2 ( x 1 y 2 − − x 2 y 1 + ⋯ ⋯ + x n y 1 − − x 1 y n ) {\displaystyle \qquad ={\frac {1}{2}}{\Big (}x_{1}y_{2}-x_{2}y_{1}+\cdots +x_{n}y_{\color {red}1}-x_{\color {red}1}y_{n}{\Big )}}

Da die Hälfte der Determinante die vorzeichenbehaftete Fläche des Dreiecks O P i P i + 1 {\displaystyle OP_{i}P_{i+1}} ist, wird diese Formel auch als Dreiecksform bezeichnet.

Setzt man P 0 = P n {\displaystyle P_{0}=P_{n}} (zusätzlich zu P n + 1 = P 1 {\displaystyle P_{n+1}=P_{1}} (siehe oben)), so gilt ∑ ∑ i = 1 n x i y i + 1 = ∑ ∑ i = 1 n x i − − 1 y i {\displaystyle \ \sum _{i=1}^{n}x_{i}y_{i+1}=\sum _{i=1}^{n}x_{i-1}y_{i}\ } und

2 A = ∑ ∑ i = 1 n ( x i y i + 1 − − x i + 1 y i ) = ∑ ∑ i = 1 n x i y i + 1 − − ∑ ∑ i = 1 n x i + 1 y i = ∑ ∑ i = 1 n x i − − 1 y i − − ∑ ∑ i = 1 n x i + 1 y i {\displaystyle 2A=\sum _{i=1}^{n}(x_{i}y_{i+1}-x_{i+1}y_{i})=\sum _{i=1}^{n}x_{i}y_{i+1}-\sum _{i=1}^{n}x_{i+1}y_{i}=\sum _{i=1}^{n}{\color {red}x_{i-1}y_{i}}-\sum _{i=1}^{n}x_{i+1}y_{i}}

Führt man beide Summen wieder zusammen und klammert y i {\displaystyle y_{i}} aus, so erhält man eine weitere Darstellung der Flächenformel:cite-ref-2[2]

A = 1 2 ∑ ∑ i = 1 n y i ( x i − − 1 − − x i + 1 ) {\displaystyle A={\frac {1}{2}}\sum _{i=1}^{n}y_{i}(x_{i-1}-x_{i+1})}

= 1 2 ( y 1 ( x n − − x 2 ) + y 2 ( x 1 − − x 3 ) + ⋯ ⋯ + y n ( x n − − 1 − − x 1 ) ) {\displaystyle \quad ={\frac {1}{2}}{\Big (}y_{1}(x_{\color {red}n}-x_{2})+y_{2}(x_{1}-x_{3})+\cdots +y_{n}(x_{n-1}-x_{\color {red}1}){\Big )}}

Verwendet man ∑ ∑ i = 1 n x i + 1 y i = ∑ ∑ i = 1 n x i y i − − 1 {\displaystyle \ \sum _{i=1}^{n}x_{i+1}y_{i}=\sum _{i=1}^{n}x_{i}y_{i-1}\ } erhält man

A = 1 2 ∑ ∑ i = 1 n x i ( y i + 1 − − y i − − 1 ) {\displaystyle A={\frac {1}{2}}\sum _{i=1}^{n}x_{i}(y_{i+1}-y_{i-1})}

Geht man von einem Polygon mit negativer Orientierung aus, ist auch der Flächeninhalt A {\displaystyle A} negativ.

Hinweis: In der Geodäsie ist die x-Achse vertikal und y-Achse horizontal und die Orientierung des Polygons entgegengesetzt. Vertauschen der Koordinaten und der Orientierung des Polygons bewirken bei der Anwendung der Formeln keine Änderung. Falls man eine dieser Änderungen nicht beachtet, erhält man mit den obigen Formeln in jedem Fall mit | A | {\displaystyle |A|} den gesuchten Flächeninhalt.

Die beiden letzten Formeln werden in Büchern über Vermessungskunde auch als Gaußsche Dreiecksformeln bezeichnet.cite-ref-3[3]cite-ref-4[4]

Speziell für polygonale Flächen mit Gitterpunkten als Ecken lässt sich der Satz von Pick anwenden. Andere Flächen lassen sich in der Regel problemlos durch Polygone approximieren, so dass man leicht an einen Näherungswert kommen kann.

Für das 5-Eck mit den Punkten

P 1 = ( 1 , 6 ) , P 2 = ( 3 , 1 ) , P 3 = ( 7 , 2 ) , {\displaystyle P_{1}=(1,6),P_{2}=(3,1),P_{3}=(7,2),}

P 4 = ( 4 , 4 ) , P 5 = ( 8 , 5 ) {\displaystyle P_{4}=(4,4),P_{5}=(8,5)}

ergibt sich

2 A = | 1 3 6 1 | + | 3 7 1 2 | + | 7 4 2 4 | + | 4 8 4 5 | + | 8 1 5 6 | {\displaystyle 2A={\begin{vmatrix}1&3\\6&1\end{vmatrix}}+{\begin{vmatrix}3&7\\1&2\end{vmatrix}}+{\begin{vmatrix}7&4\\2&4\end{vmatrix}}+{\begin{vmatrix}4&8\\4&5\end{vmatrix}}+{\begin{vmatrix}8&1\\5&6\end{vmatrix}}}

= 1 − − 18 + 6 − − 7 + 28 − − 8 + 20 − − 32 + 48 − − 5 = 33 {\displaystyle \quad \ =1-18\;+6-7\;+28-8\;+20-32\;+48-5=33}

A = 16 , 5 {\displaystyle \ \;A=16{,}5}

In der englischen Literatur gibt es ein Schema, das das Berechnen der 2x2-Determinanten optimiert: Das Schnürsenkel-Schema (engl. shoelace formula) (siehe Bild). Diese plastische Beschreibung zeigt die praktische Bedeutung der Gaußschen Trapezformel. Statt 10 Spalten genügen bei dieser Methode 6 Spalten.

cite-note-11. ↑ Max Koecher, Aloys Krieg: Ebene Geometrie. Springer-Verlag, 2013, ISBN 978-3-662-06809-0, S. 116.

cite-note-22. ↑ Josef Schlesinger: Der Tachygraph. Centralblatt für das gesamte Forstwesen: Organ der K.K. Forstlichen Versuchsanstalt in Mariabrunn, Band 2, Wien, 1876, S. 243.

cite-note-33. ↑ Martin Näbauer: Vermessungskunde. Springer-Verlag, 2013, ISBN 978-3-662-41866-6, S. 341.

cite-note-44. ↑ Heribert Kahmen: Vermessungskunde. Walter de Gruyter, 2020, ISBN 978-3-11-087406-8, S. 259.

• P. Bender, Uni Paderborn: Eine einfache Formel für den Flächeninhalt von Polygonen

• Beat Brüderlin, Andreas Meier: Computergrafik und Geometrisches Modellieren. Springer-Verlag, 2013, ISBN 978-3-322-80111-1, S. 36.

• P. Grobstich, G. Strey: Mathematik für Bauingenieure. Springer-Verlag, 2013, ISBN 978-3-322-80051-0, S. 113 (Dreiecksregel)

• Pietro Labranca: Probleme der Festigkeitslehre: Berechnung der Querschnittswerte und der Spannungen. Springer-Verlag, 2013, ISBN 978-3-663-13976-8, S. 69.

• Kleine Enzyklopädie Mathematik. Harri Deutsch Verlag, Frankfurt, 1977, S. 318.